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Abstract 

We study analytically and numerically a model of random sequential 
adsorption (RSA) of segments on a line, subject to some constraints sug- 
gested by two kinds of physical situations: (i) deposition of dimers on 
a lattice where the sites have a spatial extension; (ii) deposition of ex- 
tended particles which must overlap one (or several) adsorbing sites on 
the substrate. Both systems involve discrete and continuous degrees of 
freedom, and, in one dimension, are equivalent to our model, which de- 
pends on one length parameter. When this parameter is varied, the model 
interpolates between a variety of known situations : monomers on a lat- 
tice, "car-parking" problem, dimers on a lattice. An analysis of the long- 
time behaviour of the coverage as a function of the parameter exhibits 
an anomalous approach to the jamming limit at the transition point 
between the fast exponential kinetics, characteristic of the lattice model, 
and the \/t law of the continuous one. 
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1 Introduction 



The model of random sequential adsorption describes deposition processes in 
which desorption is negligible and surface diffusion very slow on the experimen- 
tal time scale. Particles land successively and randomly on the surface; due to 
the geometrical exclusion effect, if an incoming particle overlaps a previously 
deposited one, it is rejected. This model applies to many physical situations 
such as adsorption of latex balls, proteins or chemisorption at low temperature 
Q. The substrate is either a lattice or a continuous surface, depending on 
the size of the particles relative to the microscopic scale. Many versions of the 
model have been studied to adapt it to various physical situations, but all share 
a common universal behaviour for the long time approach to the asymptotic 
coverage, depending on the discrete or continuous nature of the substrate. As 
long as the minimum interval between two neighbouring particles on the sub- 
strate remains non-zero, which is the case on a lattice, one can show that the 
jamming limit is approached exponentially [Q, ^. Conversely, on a continuum 
substrate, the kinetics follows a power-law decay whose exponent depends on 
the number of degrees of freedom per particle [3- 11]. For instance, in the one 
dimensional case, where the model is exactly soluble, we have 

ek{t) = ek{oo)-A{k)e-^' 

for the deposition of A;-mers on a lattice with flux (p, whereas 

e{t) = e{oo)- A/t 



for the continuous "car-parking" problem] 12 1. 



In fact, at the mesoscopic scale, many physical situations involve both con- 
tinuous and discrete degrees of freedom. For example, with the recent advances 
in nano-technologies, it is conceivable to realize a lattice of small gold droplets 
deposited on a silicon surface |13]. The adsorption on this substrate of elon- 



gated particles (polymers) whose length is of the order of the distance between 
two gold dots can be modelized by the adsorption of dimers |Q. However, since 
the size of the gold dots (typically some nanometers of radius) is not negligible 
with respect to the length of the particles, there will be a continuous range of 
positions for the dimer to be fixed on two neighbouring dots. 

Alternatively, large particles like proteins or enzymes can be absorbed on a 
latticized substrate. Such a physical situation is described in ref.[|l^] and mod- 
elized in the following way : adsorbing sites are randomly or regularly disposed 
on a continuum substrate. Extended particles, represented by disks, land on the 
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surface but remain stuck only if they overlap one (or several) adsorbing site(s). 
Here again, the deposition process is driven by the discrete degree of freedom 
imposed by the location of the sites, and by a contimioTis one associated with 
the position of the adsorbed particle with respect to the site. 

In this paper, we investigate a one-dimensional RSA model which involves 
both discrete and continuous degrees of freedom. We show that the two above- 
mentioned physical situations, when reduced to a one-dimensional substrate, 
are equivalent to our model. In section 2, we present the model and analyse its 
jamming limit using the result of a numerical simulation. Section 3 is devoted 
to the analysis of the long time behaviour of the coverage as derived from the 
master equations of the model. Our results and the conclusions are summarized 
in the section 4. The Appendices contain the technical details of the derivation 
of the master equations (Appendix A) and of their solution (Appendix B). 

2 The model 

2.1 Two equivalent one-dimensional models 

Consider the first example depicted in the introduction. It can be considered 
as a generalisation of the lattice dimer model, in which the lattice sites have a 
non zero extension which we set to unity. In the deposition process, the dimer 
ends must stick on this discrete set of continuous intervals, the situation where 
a given interval contains both the end of one dimer and the origin of another 
being now allowed (see figure 1(a)). Moreover, the dimer length is no longer 
constrained to be equal to the lattice spacing, but is arbitrary within bounds 
compatible with these adsorbing rule, thereby introducing a new length scale. 
We denote by a the edge-to-edge distance between two neighbouring sites. The 
dependence of the model on this parameter occurs only in a trivial redefinition 
of the effective flux of incoming particles. Therefore, the model is independent 
of a, provided we impose the condition that a deposited particle overlaps an 
inter-site interval. The length i of the dimer lies in the interval [a, a + 2]. We 
set £ = a + r, where r G [0, 2] will be the unique parameter of the model. In 
the following, we will refer to this model as model (I). 

Consider now the second example of the introduction where the adsorbing 
sites are regularly disposed and the radius of the particles is smaller than the 
inter-site distance. In one dimension, the substrate is a regular lattice, the 
particle a segment of length r < 2 lattice units, and the adsorbing rule, aside 
from the RSA ones, is that the segment must cover a site (see figure 1(b)). We 
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call this model, model (II). 

If we constrain the deposited particle to overlap only one adsorbing site, it is 
obvious that model (II) is equivalent to the previous one taken in the limit a = 0, 
where the inter-site interval, reduced to a point, is identified with the adsorbing 
site of model (II). In the following, we will consider the one dimensional model 
from this point of view and from now on everything will refer to model (II). 

2.2 The jamming limit 

The constraint that a deposited particle must overlap one and only one adsorb- 
ing site just modifies the incoming flux (/> by a multiplicative factor of r for 
r < 1 and 2 — r for r > 1. We define the occupancy rate as the fraction of 
occupied sites. The covering of the whole substrate is rO. 



The jamming limit for this quantity, obtained from a numerical simulation, 
is presented in figure 2, as a function of r. Let us first analyse this curve qualita- 
tively. Clearly, for r < 1/2, the adsorption of a particle on the site i, regardless 
of its position on the site, cannot prevent the deposition on neighbouring sites. 
Therefore, each site will be occupied independently of its neighbours, and the 
model is completely equivalent to the deposition of monomers on a lattice; we 
expect the asymptotic limit = 1 to be reached exponentially fast. 

Consider now for the site i the most defavourable situation in the case r < 1, 
(fig. 3(a)), where the site i — 1 is occupied by the extreme left edge of a particle 
and the site i -|- 1 by the extreme right edge of another one. In this situation 
the interval for a deposition at site i is minimum and has extension 2(1 — r). 
The deposition is allowed if 2(1 — r) > r <^=^ r < 2/3. Therefore, for r < 2/3, 
we still expect a coverage of all the sites, which is effectively observed in figure 
2 where 9{t = oo) = 1 up to r = 2/3. Furthermore, since in this r-interval the 
smallest target is of non-zero extension, the kinetics must remain lattice-like 
and the jamming limit is again approached exponentially fast. 

Conversely consider now the most favourable situation in the case r > 1 

(fig. 3(b)), where the particle deposited on the site i — 1 has its extreme left 

edge close to site i — 2. The remaining space to adsorb a particle on site i, 

has extension 3 — r. If this interval is less than r (=^ r > 3/2), deposition is 

impossible, and the particle at site i — 1 effectively occupies two sites Q In this 

situation, for r > 3/2, the model is completely equivalent to the lattice dimer 

1 - e-2 

model [|16[ and one expects 9{t = oo) = = 0.43233... (see figure 2). 



^This result is more clearly seen in the framework of model (I) 
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Between these two extreme cases, for 2/3 < r < 3/2, the asymptotic oc- 
cupancy rate 9{t = oo) decreases continuously. Furthermore, since the target 
intervals for a particle deposition can be arbitrarily small, one expects a power 
law dynamics, characteristic of a continuous model. The case r = 1 is special: 
the particle length matches exactly the inter-site distance and the constraint for 
a landing particle to overlap one and only one adsorbing site is automatically 
satisfied; the discrete nature of the substrate no longer plays any role and we 
recover the "car-parking" problem|l^ with jamming limit 6 = 0.747597.... 



3 The long time behaviour of the coverage 
3.1 The master equations 

We have derived the master equations for the time evolution of the probability 
of finding at time t a gap of given length. All the technical details are given in 
appendix A. The result is that the coverage 9{t) can be expressed in terms of a 
reduced probability P{x,t) : 

e{t) = f dt' e"^' r dx P{x,t') P{r - x,t') for r g]0, 1] (1) 
Jo Jo 

= f dt' e^^-''^^' f dx P{x,t') P{r-x,t') forrG[l,2[ (2) 

Jo Jr-l 

where P{x, t) is solution of the integral equation 
case r < 1 x € [0, rl 



xP{x,t)+ r P{x',t)dx +u{x + r -I) e'"'^ T^"^ ^ P{x' ,t) 

Jx Jo 



_dP{x,t) 
dt 

case r > 1 x€[r — 1,1] 



dt 



(x-r + l)P(2;,t) + /^P(x',t)(ix' + e-(2-'-)* t P[x' ,t)di^\) 

Jx Jr-l 



subject to the initial condition P(x,0) = 1. In eq.(^), u{t) is the Heaviside step 
function and in the last integral of eq.(Q), h[x) = Inf {x -|- r — 1, 1}. 

These equations can be exactly solved in a few special cases, corresponding 
to r < 1/2, r = 1 and r > 3/2, which are presented in the next section. For the 
generic case 1/2 < r < 3/2 (r ^ 1), we have devised an iterative construction 
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of the general solution, which is developed in appendix B and which leads to 
the following result : 

Due to the last term in eqs(^,§), P{x, t) appears to be piecewise defined with 
respect to x, in successive intervals of length (1 — r) within the range (0, r) for 
r < 1, or of length (r — 1) within the range (r — 1, 2 — r) for r > 1. The number 
of intervals thus depends on the value of r and is given by (see appendix A) : 



k 



1 



+ 1 forl/2<r<l 



k 



1 



+ 1 for 1< r < 3/2 (5) 



where [X] denotes the integer part of X. 

The function Pi{x,t) which coincides with P{x,t) in the interval is ex- 
pressed in terms of a unique x-independent function q{t) 

Pf,{x,t) = q{t)+ t Ki,{x,t\t') q{t') dt 2<£<k (6) 
Jo 

Pl{x,t) = q{t) 

where the kernels K£{x,t\t') are constructed recursively. The function q{t) is 
itself solution of a linear integral equation 

q{t)=Uo{t)S^l + Pkit,t') q{t') dt'l (7) 

where mo(^) is defined by 

l/2<r<l uo{t) = e-''^ 
1 <r < 3/2 uo{t) = exp 



<2~r)t _ (2 _ r)t - 1 



and the kernels pk are expressed in terms of K£. 



Even in the simplest case, k = 2, we are unable to find analytic solutions of 
equation @. However, from these expressions we can extract both approximate 
solutions and the exact asymptotic behaviour for the dynamics of the model, 
which allow us to understand the various regimes which interpolate between 
the special cases depicted in section 3.2. 
This is the object of the following sections. 

3.2 Special cases 

There are 3 special cases where the rate equations reduce to a simple form 
which is exactly solvable. 

• 0<r<l/2. The number k of intervals is equal to one (see eq.(^)). 
Equation (^) yields directly 

P{x,t) = e""* 
e{t) = l-e-'' 
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which is, as expected, the exact result for the deposition of monomers on 
a lattice, up to a factor r in the flux corresponding to the target area for 
each deposition. 

• r = 1. The two rate equations (|3|,^) become identical and reduce to the 
known car-parking equation under the ansatz P{x,t) = e~^* Q{t)- 

• 3/2<r<2. There is again only one interval. From the rate equation we 
get directly 

P{x, t) = uo{t) = exp [e~(^-'^)* - (2 - r)t - 1 

and 

9{t) = (2 - r) /* dt' e^^-'-)*' ul{t') = - [l - exp[-2 + 2e-(2--)*] 
Jo 2 L 

which is exactly the lattice dimer solution, up to a factor 1/2 due to our 
definition of 9. 

3.3 The case 1/2 < r < 2/3 

According to eq.(|5[), this is the interval (for r < 1) corresponding to the lowest 
non trivial value k = 2. There are only two reduced probabilities, q{t) and 
P2{x,t). The function p2{t,t') of eq.(^) is given by 

P^(t^*') = - (2^-1) 

and K2{x,t\t') of eq.(^) by 

K2{x,t\t') = e-'' ^-^^ 

(see the appendix). After some elementary algebra based on the integral rep- 
resentation eq.(|6[), one gets from eq.(|l|) a very simple expression for 9 : 

9 = l-q\t)e''' (8) 

We obtain the exact asymptotic behaviour of the coverage in the following way. 

From the inequality < q{t) < 1 and equation it follows that q{t) < C ^ 

which implies that the functions 

t 



Cn{x,t) = / e-("-^)*' q{t') t'" dt' 
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have a positive, finite, t ^ oo limit for all n > and for all x € [0,r], that we 
denote Cn{x). Expanding the denominator of the kernel p2{t,t') in power of 
t' /t, we get from eq.(0) the asymptotic expansion of q{t) at large t 



lit) 



-{l-r)t 



C„{l-r) + ^l(l^ + 0(^; 



(9) 



Inserting this result in eq.(P), we obtain for the coverage the asymptotic value 
6{t = oo) = 1, as expected, with the approach : 

-(2-3r)i 



9{t) = l-C^{l-r) 



t2 



+ 



(10) 



This equation shows that the kinetics remains exponentially driven over the 
open interval 1/2 < r < 2/3 but becomes l/t^ at the end point r = 2/3. 

3.4 The case 2/3 < r < 3/4 

We have explicitly studied this case, where three x-intervals {k = 3) are in- 
volved, since we expect it to be typical of what happens over the remaining 
range r < 1. 



Following the method of the previous section, we can derive for the reduced 
probabilities P2 and P3 the following large-time behaviour 



P2ix,t) = Coix,t) 



-xt 



t 



P^{x,t) = P2{x,t)+0 



l + 0(^) 

g-(2r-l)A 



whereas the asymptotic behaviour of q{t) is still given by eq.(|9|) in spite of its 
being defined from a different kernel p^. 

The time derivative of 6 may be expressed in terms of the reduced probabilities 
as follows 



e{t) = e 



rt 



2r-l 



1-r 



P2{x,t) P2{r-x,t) dx + 2q{t) 



2{l~r) 



2r-l 



P2{x,t) dx + 2q{t) 



2(l-r) 



^3(2;, t) dx 



Inserting in this expression the previous asymptotic expansion, we get 
where A and B are constants depending only on r. In particular 

f.2r-l 



(11) 



A = lim 



1-r 



Co(x, t) Co{r — X, t) dx 
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We see that, except at r = 2/3 where A vanishes, the leading behaviour of 9 is 
of the form 

e{t) = 0(oo) - ^ 

We expect this behaviour to hold over the whole range up to r = 1 where it is 
proved. We have checked this assumption numerically for r = 0.70 and r = 0.9. 
For r = 2/3 we recover the result of the previous section 

9{t) = ^(oo) - I 

which confirms that this unusual behaviour occurs only at this special value of 
r. 



p,{t,t') = e-(2-)*' / ' + L " ^ ^ " exp [l - e~(^-)(*'+^")j dt" 



3.5 The case 4/3 < r < 3/2 

As in the previous section, this interval corresponds to the first non-trivial one 
where A; = 2 in the case r > 1. The kernels K2 and p2 are defined by 

K2{x,t\t') = e-*'e('-i)*^ 

rt-t' e-(3-2r)q^(3 _ 2r)t"] - 1 
/o t 

The expression for 6 in terms of the reduced probability q{t) is obtained from 
eq.® : 

6 = 1- q\t) e2(2--)* - (2 - r) /* e^^"'^)*' q\t') dt' (12) 

JO 

One observes that p2 is negative, which implies the bound < q{t) < uo{t) and, 
from the mean-value theorem, the estimate 

P2(M') = G{t-t') I ' (3-2r)| 

where 1 < exp(l - e'^^-^)*') < G{t - t') < exp(l - e-^^-*^)*) < e = 2.718.... In- 
serting this result in eq.(|^ leads to the asymptotic time expansion of q{t) 

qit) = e-(2-'-)* \A{r) + ^ + 



where the bracket reduces to a constant when r goes to 3/2. This shows from 
eq.(p!^) that the jamming limit 9{oo) is reached as 1/t. 
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4 Summary and conclusions 



We have studied a one- dimensional model of random sequential adsorption with 
discrete and continuous degrees of freedom, which depends on one length pa- 
rameter r. When this parameter varies over its range < r < 2 the model goes 
through the following regimes : 

• < r < 1/2 

monomer on a lattice; 9{t) = 1 — e"*"* 

• 1/2 < r < 2/3 

total coverage; non trivial lattice dynamics : 6{t) = 1 — Ae~^'^~^'^^^ /t'^ for 
large t. 

• r = 2/3 

total coverage; anomalous continuous dynamics : 9{t) = 1 — A/t"^ ioi large 
t. 

• 2/3 < r < 3/4 

non-trivial asymptotic coverage; normal continuous dynamics : 9{t) = 
6{oo) - A/t for large t. 

• r = 1 

"car-parking" problem : 9{t) = 0.7476 ■ ■ ■ — A/t for large t. 

• 4/3 < r < 3/2 

same as for 2/3 < r < 3/4. 

• 3/2 < r < 2 

lattice dimer model; 9(t) = ^ 1 - exp[-2 + 2e"(^"'")*] 

• 3/4 < r < 4/3 

Although we have not investigated this interval analytically (except for 
r = 1), we expect the same regime as in the bordering intervals 2/3 < 
r < 3/4 and 4/3 < r < 3/2 to hold, that is a non trivial asymptotic 
coverage decreasing with r and a 1/t normal continuous dynamics. We 
have checked this conjecture numerically. 

The remarkable point is that the kinetics of the model exhibits three "phases" : 
for < r < 2/3 where it is lattice-like, for2/3<r<3/2 where it is continuous- 
like and for 3/2 < r < 2 where it is lattice-like again. At the transition point 
r = 2/3 it becomes "anomalous" since the jamming limit is approached in l/t^ 
in contrast with the general belief that the exponent n of the power-law decay 
is equal to the inverse of the number of degrees of freedom per particle. 
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The regime on both sides of the transition is characterised by a same typical 
cross-over time 1 17] r = \ 2-3r\ ' defined from the slope of the exponential in 
eq.(|lO|) or from the ratio B/A in eq.(ll). This time is such that for t <^t the 
dynamics is dominated by a 1 /t^ behaviour in the two "phases" ; only for t ^ t 
does the characteristic long time behaviour, exponential for r < 2/3, in 1/t for 
r > 2/3, emerge. Since r — > oo when r — > 2/3, this long time regime is squeezed 
at the transition point r = 2/3, leaving only the 1/t^ behaviour. 

To find for 9{t) an expression in a closed form is a very difficult techni- 
cal problem, comparable to the determination of the correlation function in 
the standard "car-parking" model|18]. However, the properties of the kernel 
yOfc allow us to obtain an iterative solution of eq.(|7|) from which approximate 
expressions for the long time coverage can be obtained. 

Finally, it may be interesting to investigate these kinds of models in a more 
realistic physical context, such as a two dimensional substrate, a disordered 
distribution of adsorbing sites or the possibility for a particle to overlap several 
sites. It is presumably difficult to get an analytical insight in such models, but 
their properties can be numerically analysed using as a guide the one dimen- 
sional results. 
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APPENDIX 



A The master equations 

In this appendix, wc derive the master equations for the probability of finding 
at time t an unoccupied interval (gap) of given length. Let Pn{x,y,t) be the 
probability for finding a gap of length at least x + y + n — 1, where n is the 
number of adsorbing sites in the gap and x{y) the distance between the left 
(right) edge of the gap and the last left (right) site in the gap (see figure 4). 
Pi(0,0,t) is the probability for finding a gap containing at least one site and 
its complement, 1 — Pi(0, 0,t), defines the probability that a site is covered by 
a particle, which is equivalent to the occupancy rate : 

e{t) = 1 - Pi{o,o,t) 

The rate equation for this quantity is (in these equations and in the subsequent 
rate equations we set the effective flux of particles to unity) 

_m(O^M ^ [' PAx,r-x,t)dx for re [0,1] (13) 
ot Jo 

= I Pi{x,r - x,t) dx forrG[l,2] (14) 

Jr-l 

Therefore to determine 9 one needs to know Pi(x,y,t) only for x and y less 
than r (case r < 1) or greater than r — 1 (case r > 1) and set x + y = r. The 
general rate equations are obtained as usual, by counting the different ways of 
destroying a gap : 

case r < 1 x G [0, r] y & [0,r] n >2 
dPn{x,y,t) _ 

at 



(15) 



[x + y + {n-2)r]Pn{x,y,t)+ f Pn{x' ,y,t)dx' + [ Pn{x,y' ,t)dy' 

Jx Jy 
rx+r—1 ry+r—l 

+u{x + r-l) Pn+i{x',y,t)dx' + u{y + r -1) Pn+i{x,y' ,t)dy' 

Jo Jo 



case r > 1 x G [r — 1, 1] y G [r — 1, 1] n > 2 

dPn{x,y,t) 



dt 



[(x - r + 1) + (y - r + 1) + (n - 2)(2 - r)]P„(x, y, t) 



+ / Pn{x',y,t)dx' + [ Pn{x,y',t)dy' 

Jx Jy 

Mx) My) 

+ / Pn+i{x',y,t)dx' + Pn+iix,y',t)dy' (16) 

Jr-l Jr-l 
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where u{t) is the Heaviside step function and, in the last equation, b{x) = 
Inf {x + r — 1, 1}. The probabiUty Pn{x, y, t) is subject to the initial condition 

P„(x,y,t = 0) = l 

From these equations, we see that for n >2, the n dependence can be factorised 
out : 

Pn{x,y,t) = e-^''-^>' P2{x,y,t) for r G [0, 1] (17) 
= e-("-2)(2-r)t p^^^^y^^) forrG[l,2] (18) 

The function P2{x,y,t) must be symmetric in its spatial arguments. Moreover, 
for n > 2 the left and the right parts of the gap cannot be simultaneously 
affected by the deposition of one dimer, hence the x and y dependence are 
uncorrelated. We deduce from these considerations that : 

P2{x,y,t) = P{x,t) Piy,t) (19) 

reducing the problem to finding one unknown function, P(^x,t), which we call 
in the following, a reduced probability. 

The rate equation for the function Pi{x,y,t) is different from eqs( |T^jT^ ) and 
depends on the value of x + y with respect to r. However, the arguments leading 
to the factorisation of the x and y dependences and of the n dependence remain 
valid for x + y > r which is precisely the region of interest, leading to : 

Pi{x,y,t) = e'' P{x,t) P{y,t) x + y>r 

From this factorisation property and the rate equations (|l3[0), we can express 
the occupancy rate in terms of the function P, leading to eqs.(|l],||) : 

e{t) = f dt' e'"*' dx P{x, t') P{r - x, t') for r G [0, 1] 
Jo Jo 

= f dt' e(2-^)*' / dx P{x, t') P{r - x, t') for r G [1, 2] 

Jo Jr-l 

Prom eqs. (|T5| , |l^ ) and the factorisation properties of P„(a;, y, t), eqs.(^, 18,19) 
we deduce the rate equations for the function P, eqs.(|^,^) of section 3.1 : 

case r < 1 x G [0, r] 



dP{x,t) 
dt 

case r > 1 x G [r — 1, 1] 



xP{x,t)+ r P{x',t)dx' + u{x + r -I) e""^ T^'' ^ P{x' ,t 

Jx Jo 



.^^^lll = {x-r + l)P{x,t)+ P{x',t)dx' + e-^^~''^' f^''''^ P{x',t)42i) 

Ot Jx Jr-l 
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Due to the last term in eqs.(^,21), P{x, t) appears to be piecewise defined with 
respect to x, in successive intervals of length (1 — r) within the range (0, r) for 
r < 1, or of length (r — 1) within the range (r — 1, 2 — r) for r > 1. To be more 
explicit, let us define the set of intervals le in the following way : 



k — 1 

For < r < 1 define k > 1 such that — : — < r < 



k 



k - k + V 

I,= [{£-l){l-r),£il-r)] 

4 = P-l)(l-r),r] 
k 

which are such that [J -^^ = [0) '''] 



then : 

l<i<k-l 



k + 2 k + 1 

For 1 < r < 2 define k > 1 such that < r < ; then : 

k + 1 ~ k 

Ii = [{2-r),l] 

onlyfor/i;>2 h = [{k - l){r - 1), {2 - r)] 

only for /c> 3 le = [{k - £ + l){r - I), {k - i + 2){r - 1)] 3<£<k 

k 

which are such that U = ~ li 1] 

1=1 

The number k of such intervals is directly expressed in terms of r : 



k 



l-r 



+ 1 for < r < 1 



k 



1 



+ 1 for 1 < r < 2 



where [X] means integer part of X. In the following, for a given value of r, we 
will denote by Pe{x, t), ^ = 1, k the restriction of P{x, t) to the interval Ig. 



B Integral equations for the generic case 2/3 < r < 

3/2, r + 1 

In this appendix, we derive the representation of eqs.(^J^ from eqs(^,21) and 
we give the expression of the kernels Kg and /O^ and of the function uq. 



We first remark that the reduced probability Pi{x,t) is in fact independent 
of X : Pi{x,t) = q{t). Consider the left bordering site of the gap, in the case 
r < 1. There is a vacant space of length at least 1 — r from its right edge 
independent on the occupation of its left neighbour; this means that P{x, t) is 
independent of x for x £ [0,1 — r] . The same property holds, in the case r > 1, 
for x £ [2 — r, 1]. 
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We start with the case r < 1. For 1 — r < x < r, we derive eq.(20) with 
respect to x to get 



dxdt dx 



+ e-''^P{x + r-l,t) 



dP(x, 0) 

Then integrating with respect to t, using the initial condition — ^ — = 0, we 

ox 

obtain 

= - /* e-*'"-(*-*') P{x + r- l,t') dt' 
ox Jo 

Finally, we integrate with respect to x, taking into account the boundary con- 
dition P{x = 1 — r,t) = q{t), to get 

P(x, t) = q{t) - e-(i^'-)* /* dt' e-(2'-i)*' dx' e"' P{x', t') (22) 

JO JQ 



By using repeatedly eq.(22) as x increases, one obtains equation (^) : 

(x, t) = q{t) + / Ki{x, t\t') q{t') dt 2<i<k 
Jo 

For X G I2 then x' G Ii where P{x',t') = q{t') and eq.(p2|) gives eq. d) with 
^ = 2, and the following expression for the kernel K2 

-x{t-t') _ -(l-r)(t-t') 

K2{x, t\t') = e-^* r < 1 (23) 



When X S Is, by including the previous result in eq.(22), we get eq.(^) for 
= 3 with the explicit form of the kernel i^3, and so on. 



Analogous manipulations can be performed on eq.(]2l]) in the case r > 1 for 
X < 2 — r, where the boundary condition is now P{x = 2 — r,t) = q(t).Th.e 
equation equivalent to eq.(p2|) is 



P(x, t) = qit) + e2('-i)* /* dt' e-"*' dx' e^''(*'-*) P(x', t') (24) 

Jo Jx+r-l 

which, by repeated application as x decreases, gives the representation of eq.(^). 
For x G I2, one gets the kernel K2 



e Mt-t') _ g-(2-r)(t-t') 



i^2(x,t|t') = e-*'+('-i)* r > 1 (25) 



Assuming that, for a given value of r all the are known, one can define 
a kernel crk{t, t') 

(^k{t,t') = ^ / dx Ke{x,t\t') 
1=2 -^^^ 
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and derive an integral equation for q{t). Considering eqs. (M P) where X IS 
fixed to the value x = 1 — r (case r<l)orx = 2 — r (case r > 1) and using 
eq.(|6|) for Pi{x,t), one gets respectively : 



r < 1 
r > 1 



^(t) = rq{t) + J\{t') ak{t,t') dt' 
dq 



dt 



(t) = {2-r)[l + e-(^-^^']q{t) + e-^^"-^' qit')ak{t,t')dt' 



By integrating with respect to t with the initial condition q{0) = 1, we get the 
equation (0) 

(?(t) =Uo(t)(l+ / Pk{t,t') q{t')dt' 



with 

r < 1 UQ{t) = e"''* 
r > 1 lio(t) = exp 



-(2-r)t 



{2-r)t-l , pfc(f,0 



, Pk{t,t')=J^ e^'" ak{t",t')dt" 



exp 



1 - e 



-(2-r)t"" 
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Figure captions 



Figure 1 (a) description and parametrization of model (I) : deposition of 
dimers on extended adsorbing sites, (b) description and parametrization 
of model (II) : deposition of segments on localised adsorbing sites. 

Figure 2 The jamming limit of the model as a function of the length r of 
the adsorbed segments. The values of r separating the different regimes 
are indicated. 

Figure 3 Two extreme situations for the deposition in model (II). 
Figure 4 The parametrisation of the gap. 
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